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In this paper, the approach for investigation of asymptotic {Re — ► oo) scaling exponents of Eulerian 
structure functions (J. Schumacher et al, New. J. of Physics 9, 89 (2007). ) is generalized to 
studies of Lagrangian structure functions in turbulence. The novel "bridging relation" based on 
the derived expression for the fluctuating, moment-order - dependent dissipation time T v , n , led to 
analytic expression for scaling exponents («„ ) of the moments of Lagrangian velocity differences 
S n ,L(f) = (u(t + t) — u(t)) n oc r Kn in a good agreement with experimental and numerical data. 



Introduction. A turbulent flow can be described us- 
ing Eulerian and Lagrangian approaches addressing dy- 
namics of velocity field and evolution of individual fluid 
particles, respectively. Therefore, one can introduce 
two kinds of structure functions, i.e. moments of ve- 
locity increments. The properties of Eulerian corre- 
lation functions (ESF) were first theoretically investi- 
gated by Kolmogorov in his celebrated K41 theory of 
turbulence based on an exact in the inertial range re- 
lation for the third-order structure function S^r) = 
(u(x + r) - u(x)) • £) 3 = (<5 r u • £) 3 cx r [1]. The prob- 
lem of Lagrangian correlation functions, for which not 
a single exact dynamic relation exists, is not new and 
was originally formulated within the framework of Kol- 
mogorov theory of turbulence (for a detailed review, see 
Ref. [2]). Due to technological limitations of the past, 
experimental studies of strong turbulence were mainly 
devoted to Eulerian structure functions (ESF) S n — 
(u(x + r) — u{x)) n = (S r u) n limited to the single- point 
measurements with subsequent application of Taylor hy- 
pothesis. Here u is a component of velocity field parallel 
to the displacement vector r chosen along the rc-axis. The 
experimental and numerical studies of Eulerian structure 
functions revealed two distinct intervals: while in the an- 
alytic range r — > (or r <C r/„), the velocity field is 
differentiable and S n {r) oc r™, at the scales L»r>?|„, 
the ESFs are given by algebraic relations S n (r) oc r^™ 
with 'anomalous' scaling exponents £„. By this defini- 
tion, ?7 n 's are the cut-offs separating analytic and 'rough' 
ranges of ESFs. (See Fig. 1). 

The remarkable breakthroughs in experimental particle 
tracking, led by the Bodenschatz group [3]- [7], enabled 
investigation of Lagrangian structure functions (LSF) 
S n ,L = (u(t + t) — u(t)) n = (6 T u) n which are a cru- 
cial ingredient in understanding of turbulent transport 
and mixing. Due to analyticity of velocity field, as 
t ► 0, S T u ~ clt where a is acceleration of a fluid 
particle and S u l — a n T n . For the time-increments 
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r r^ n, the LSFs S n i oc r K ™ with the dissipation 
times Tv^n separating analytic and rough intervals on 
a time-domain. (See Fig.l). The recent multifractal 
theory of LSFs [6] stressed the importance of the dis- 
sipation time t v which was treated in the spirit of Kol- 
mogorov theory as a moment-number- independent quan- 
tity oc T j\J Re — const where T w L/u rms is a large- 
scale eddy turn-over time. To make a connection between 
ESFs and LSFs, the authors of Ref. [6] used the "bridg- 
ing relation" (BR): r sa t5 t u = r(u(<+r) — u(t)) with the 
time-increment r approximately equal to the fluctuating 
eddy turn-over time in the inertial range. This expres- 
sion, introduced on dimensional grounds in Ref. [8], can 
be understood as follows. Consider a fluid particle at a 
time t = to occupying position Xo. Then, the particle 
displacement is: R = X(io + r) — Xo = f t ° T u(A)dA, so 

that R = u(t )T+ j* o+T (u(\)-u(t ))d\. We can see that 
if u(t) = U = const, the displacement R^O. Then, fol- 
lowing [2] (page 342, where it is related to transition to a 
moving frame of reference) , we define the regularized (not 
involving single-point- velocity) quantity: r = R— u(to)r. 
(It is estimated in Ref. [2] (page 359) that r ss 0, mean- 
ing that it tends to zero when either u —> U = const 
or r — ► 0.) From the mean value theorem we also have 
f t ° T u(X)dX w u(t + t)t, where < r < r. The relation 
r w tS t u is obtained in the first approximation setting 
u(t + t) w u(t + r) which resembles the one used in 
construction of Kraichnan's Lagrangian History Direct 
Interaction Approximation [9]. Keeping in mind topo- 
logical complexity of developed turbulence, we conclude 
that the mean value theorem, leading to the "bridging re- 
lation" r « tS t u with the inertial range time-increments 
t ^> T, q , cannot be accurate. Recently, the BR was an- 
alyzed using the exact relations between Eulerian and 
Lagrangian structure functions by Kamps et al [10] who 
showed that, combined with the multifractal formalism, 
it leads to the Lagrangian exponents in a substantial dis- 
agreement with experimental [3]-[7] and numerical [11] 
data. Moreover, it was pointed out that the theory of 
Ref. [6], expressing anomalies of Lagrangian exponents 
in terms of anomalies of Eulerian ones, does not explain 
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the " 2D-paradox" : while the two-dimensional Eulerian 
turbulence is not intermittent, the Lagrangian one is. 
In this paper, based on the ideas developed in Refs.[12]- 
[16], we attack the problem differently. In the analytic 
(dissipation) range where S T u « or, the displacement 
r w tS t u/2. Extrapolating this to the dissipation cut-off 
77, separating analytic and rough intervals of the structure 
functions, leads to the dissipation time t v : 

V~T, q S Tv U (1) 

Both rj w and t v in formula (1) are random functions 
investigated in great detail theoretically and numerically 
[13]- [16]. Unlike the BR defined for the large inertial 
range values of the time - increment t » t v , the expres- 
sion (1), introduced as an extrapolation of the exact in 
the analytic interval relation, is accurate for short times 
t < t v . Theoretical basis for this extrapolation is under- 
stood as follows. The structure functions, both Eulerian 
and Lagrangian, can be formally represented as ipn( x ) 
x Pn(x) yjjtfo th e x-dependent exponents £„ < (3 n (x) < n 
(or K n < P n (x) < n) covering both analytic (f3 n = n) and 
inertial ( (3 n {x) = or (3 n = K n ) ranges. The smallest 
inertial range exponents giving maximum values to the 
differences n — (3 n , correspond to the strongest singular- 
ities of velocity field and defining the cut - offs by the 
matching relation a n r Vtn = K n r^ n , we account for the 
strongest, dominant, singularities. (See Ref.[16]). 
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FIG. 1: Definition of the cut-off r\ n ( or r Vn ) for a structure 
function S„(x) (or(S n ,L(x)). As an example S — Ss is chosen. 
The strait lines are x 8 (analytic range) and x 2 ' 2 (inertial one). 
The point of their crossing is x n = r) n - Schematic. 

The probability distribution functions Q(n,Re), derived 
and investigated numerically in [14]- [15] enable one to 
evaluate the moments of the dissipation scales: 



rf l n = n n = J T] n Q{V, Re)dri « LReJ 



(2) 



Combining (1) and (2) and assuming continuity of S n L(r) 
on the dissipation time scale yields: 



L"ifc«»-e„ + i-i = {u rms T) n {^) 



V,n \K n +n 
T 



(3) 



Our goal now is to express the dissipation times T, q n = r„ 
in terms of the Reynolds number and, comparing the re- 
sult with formula (2), obtain expressions for Lagrangian 
exponents K n . 

Eulerian structure functions. Since energy dissipation in 
the inertial range dynamics ( inverse-energy-cascade ) of 
2D turbulence is irrelevant, the relation (2), obtained by 
balancing viscous dissipation and inertial-range contri- 
butions to the exact equations for the Eulerian struc- 
ture functions (see Refs. [14]-[16]) is valid for three- 
dimensional flows only. If, in accord with Kolmogorov 
theory, we assume = n/3, this formula gives the well- 
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known n- independent relation r\ n « LRe~ J . However, 
due to intcrmittency, the function is a convex function 
of the moment order n and for a fixed Reynolds num- 
ber Re = const, the dissipation scales r\ n decrease with 
n. (This result has been numerically tested in Ref.[16]). 
Since the inertial range is compressed to the interval be- 
tween integral and dissipation scales (r) n <C r <C L/10), 
experimental determination of exponents is very diffi- 
cult and at the present time only the exponents with 
n < 8 have been accurately established by direct inves- 
tigation of inertial range dynamics. In accord with the 
recently developed theory of small-scale intermittency, 
the dissipation scale is defined by a dynamic Reynolds 
number Re n = rj 1 [15]. First, we see that the dis- 
sipation scale is not a constant as in K41, but a fluctuat- 
ing property of a flow. Then, it is easy to show [14[-[16] 
that |j = ^ Rj iia^L. Since r\ is the scale separating 
analytic and rough scale- intervals of the velocity field, it 
has been shown that: 



(g) n » S 2n { mn )/v n ex nfr/v" cx Re"- (4) 



Using (2) we derive readily: 



Pn = n- 
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-; d n = n- 



(5) 



where d n are the exponents of the moments of the dissipa- 
tion rate £ n rj ^""^ = Si n {r)i n )/v cx Re d ™ . According 
to this calculus, the moments of Lagrangian acceleration 
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a* = (^ L -) n oc ij e «3n-«3n+i- 1+ " = R e a " (6) 

The quantitative results of Ref. [16], obtained with the 
help of a particular parametrization (Ref. [13]- [15]) 



£„ = 0.383n/(l + n/20) 



(7) 
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The predicted in Refs. [13]- [16] relations (4), (5) have 
been confirmed in the most detailed numerical simulation 
of Ref.[16] and [17]. The expression (7), calibrated to give 
£3 = 1, is the result of a theory valid for the even - order 
moments only. We do not have any reason to believe 
that \S r u\ 3 oc r^ abs with £ a b s = 1. Thus, the accuracy of 
the relation (7) must be of the order |£j — £3, &s|, which 
is numerically small. This drawback is common to all 
exisitng models leading to expression for . 
Lagrangian structure functions. The classical treatment 
of the problem can be formulated as follows [2]. As t — > 
0, the velocity field is differentiable and using (6): 



S 2 n,L(r) - (5 T u) n = a 2 "r 2 " oc Re a ^r 2n (8) 

The limit t — > 0, means that there exist a time-scale 
T~n,2n such that the relation (8) is valid in the interval 
t <C T V 2n- In the "inertial range r v , 2n -C r <C T = 
^w— ft° +T u(0)u(t)dt, the velocity field is not differen- 
tiable and 

Vi(r)ocr K - (9) 

In what follows we define n th dissipation time r 2n by the 
matching relation: 

*w(t 2 „) = a^r% = {ET 2n n T -fr n - n (io) 

According to K41, neither integral scale nor viscosity in- 
fluence the dynamics of inertial range and, on dimen- 
sional grounds n 2n — n and S 2n ,L = K 2n {£T) n with 
K n = const. Let us examine consistency of this result 
with some other predictions of the K41. It is clear 
that: S 2 ,l = fo Jq a(Xi)a(X 2 )d\idX 2 . By the time- 
homogeneity a(Xi)a(X 2 ) = A(Xi — X 2 ) = A(s) and 

S 2 ,l= I (T-s)A(s)ds (11) 
Jo 

According to K41 (see Rcf.[2]), in the interval < r < 
t v ,2, Mr) w a 2 and for t V:2 < t < T, A(t) w K 2 £/t. 
Substituting this into (11) gives the inertial range expres- 
sion (r > T Vt2 ) 

S 2 , L {T)~ a ^ + K 2 £{T\n--T + T 2 ) (12) 

2 T 2 

We can notice a slight inconsistency of K41: the expres- 
sions (12) and the K41 result S 2i l oc t differ by a large 
logarithmic factor which is impossible to derive from di- 
mensional considerations. Dimensional considerations [2] 

also give a 2 w Si/^/v oc v Re and up to logarith- 
mic correction, the second contribution to the right side 



H £n Q^n ^"n ^ In 

2 0.7 0.55 0.95 0.9-1. [4]'[10[,[11] 0.545 

4 1.27 1.5 1.44 1.3-1.6 [2], [4], [5], [11] 0.566 

6 1.77 2.6 1.74 1.6-1.8 [2], [4], [5], [11] 0.61 

8 2.2 3.9 1.92 1.9,[4] 0.66 

10 2.55 5.8 2.05 2.1, [4] 0.76 



of (12) is O(l). Demanding continuity of the structure 
function S n ,L in the limit Re — > 00, we obtain the familiar 
K41 estimate for the relaxation time t 2 sa Urms oc Re~ i . 

terms 

Accounting for intermittency. Since, T ss L/u rms , then, 
as follows from (10), on the dissipation time-scale: 

S 2 n{T 2n ) - -Re { — ) - K 2n {— ) 

tl rms 1 1 

(13) 

and for t <C r 2n : 

S 2n ,L oc i?e Q2 "T 2 "; ^ w Re^^ = Re^ (14) 

Now, we derive the expression for n n . Combining (2), (3) 
with (14) gives: 

L n Re(n -C +1 -i = (u rms T) n (^r-+ n = 

{Tu rms ) n Re a ^T-i n) (15) 

and 

n _ OnjKn + n) ^ 

£,n ~ £,n+l ~ 1 K n — n 

With the parametrization (7) (or any other, for that 
matter) one is able to calculate the exponents n n of La- 
grangian structure functions. The theoretical predictions 
(k^) are favorably compared with experimental and nu- 
merical (k^ p ) results in the Table. 

The moment-number dependence of Lagrangian and 
Eulerian exponents is presented on Fig. 2 (top curve) to- 
gether with comparison of theoretical predictions with 
experimental data of Ref. [4] . 

The non-trivial Reynolds number dependence of the re- 
laxation times t^.„, given by expression (14), is also of im- 
portance. Plotting the experimental data on S n ,l(t) in 
coordinates t/t k w T/^/Re instead of t/t„, the authors 
of Ref. [7] failed to collapse the experimental graphs for 
Sul{t, Re), which indicated a dynamic inconsistency of 
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FIG. 2: Top curve: Lagrangian and Eulerian exponents k„ 
and £ n vs. n. Lagrangian ones are "more intermittent". Bot- 
tom: comparison of predicted (formula (14)) and experimen- 
tally observed exponents k p = £ p [4]. Black line is exponents 
based on K41. 

Kolmogoorv's dissipation time. Similar phenomenon was 
described in Schumacher et. al. (Ref. [16]) dealing with 
the moments of increments Eulerian velocity field. 

Conclusions. In this paper we used the theoretically 
predicted and experimentally verified fact [16] that the 
inertial-range asymptotic exponents of Eulerian structure 
functions are closely related to the anomalous Reynolds 
-number- scaling of the moments of velocity derivatives 
defined on the fluctuating ultra-violet cut-offs given by 
expression (2). Following [15], [16], we introduced the 
moment-number-dependent dissipation time which, com- 
bined with the new "bridging relation" r\ ss t^St u, en- 
abled us to evaluate the scaling exponents of Lagrangian 
structure functions in a good agreement with experimen- 
tal and numerical data. Unlike the previously used BR, 
defined for the inertial range increments r 3> t v , this rela- 
tion, which is an extrapolation of an exact in the analytic 
range dependence r = tS t u/2, is much more accurate. 
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